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Abstract. The critical two-terminal conductance gc and the spatial fluctuations 
of critical eigenstates are investigated for a disordered two dimensional model of 
non-interacting electrons subject to spin-orbit scattering (Ando model). For square 
samples, we verify numerically the relation CTc — l/[27r(2 — -D(l))] e^//i between critical 
conductivity (Tc = Sc = (1.42 ± 0.005) e'^/Zi and the fractal information dimension 
of the electron wave function, D{1) = 1.889 ± 0.001. Through a detailed numerical 
scaling analysis of the two-terminal conductance we also estimate the critical exponent 
ly — 2.80 ± 0.04 that governs the quantum phase transition. 
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1. Introduction 

Quantum phase transitions in two dimensional (2d) disordered systems have been 
the subject of continuing research for many decades. In particular, systems of non- 
interacting particles that are invariant under time reversal, but lack the symmetry of 
spin rotation are of special importance due to their property of a complete Anderson 
transition 13 El IH El El • Thus, 2d systems with spin-orbit interaction, which belong 
to the symplectic symmetry class [2], exhibit an energy spectrum with regions of 
either localized or extended states separated by quantum critical points at which the 
correlation length diverges according to ^{E) ~ \E — Ec\~'^ with a universal exponent u. 
The energetical position of these points Ec depend on disorder strength and spin-orbit 
interaction. The corresponding critical eigenstates were found to be multi-fractal objects 
showing strong amplitude fluctuations as well as long-range spatial and energetical 
correlations [HIEIIIDI- 

Our knowledge about the critical properties like critical energy and disorder, 
or about the exponents that govern their scaling, originate mainly from numerical 
investigations of quantities like localization length 13 IH El USl CHI El or spectral 
correlations of eigenvalues which both are not directly accessible to experimental 
detection. The aim of the present paper is to present scaling results of the two-terminal 
conductance [in], a physical quantity that can easily be measured in experiment, and 
to examine a recently proposed relation between the average critical conductance and 
the generalized fractal dimensions of critical eigenstates. Therefore, we first present our 
numerical data for the two-terminal conductance and obtain the critical parameters of 
the model by a finite-size scaling analysis. In the metallic regime, the mean conductance 
is expected to increase with increasing system size while a decrease is usually found in the 
localized regime. At the critical point, the mean conductance gc should be a universal, 
system size independent quantity and equal to the critical conductivity a. Hence, the 
size scaling of our numerical data enables us to estimate the critical disorder Wc, and to 
calculate the critical conductance gc- We also find the critical exponent u which governs 
the energy and disorder dependence of the conductance near the critical point. 

Then we discuss the fractal properties of critical wave functions and calculate 
fractal dimensions D{q). Non- uniform spatial amplitude distributions of electron 
wave functions that exhibit multi-fractal properties, are a characteristic feature of the 
critical regime at metal-insulator transitions. It is well known that the corresponding 
generalized multi-fractal dimensions influence the time dependent transport properties 
[TT] . Here, we show that they also determine the conductance at a quantum critical 
point. Our main result is a confirmation of the relation D{q) = 2 — q/[j347r'^gDh] between 
the fractal dimensions and the conductivity (per spin direction), a = e^gD, proposed 
by Fal'ko and Efetov ^2]. We also examine the validity of Janssen's formula ITTIj 
which relates the Lipschitz-Holder exponent of the maximum of the multi-fractal 
distribution /(a) to the scaling parameter Ac = l/[7r(a;o — d)] for quasi-ld systems. We 
find, however, our result for the symplectic Ando model {d = 2) to be at variance with 
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this proposal. 

2. The model and methods 

2.1. Hamiltonian 

We study the two dimensional (2d) Ando model ^ defined on a square lattice by the 
Hamiltonian 

r[ = wY^ e„4c„ + ^ Kn'4c„,. (1) 

n [nn'] 

The En are uncorrelated random on-site energies chosen according to a box probability 
distribution in the range < 1/2. The parameter W measures the strength of the 
disorder. The spin dependent hopping terms between nearest neighbour sites in the 
direction parallel and perpendicular to the current fiow are 

M| = ^( ' 'V and \4 = \/| (2) 
\ —s c 1 \ —IS c I 

with + = 1. The parameter s determines hopping accompanied with a change of 
the spin of the electron. As usual, it is set to s = 0.5. The special case s = corresponds 
to a system with orthogonal symmetry which exhibits no metal- insulator transition in 
2d. Energies are measured in units V = 1 and lengths in units of the lattice constant 
a = 1. 

It is well known that for Fermi energy Ep ^ 0, the Ando model exhibits a metal- 
insulator transition as a function of disorder aX W = Wc ~ 5.8 jHl E]- The critical 
exponent for various 2d systems with symplectic symmetry has been estimated in 
numerous works using different numerical methods O [13 1201 HH 1221 Ull 12^1 , leading 
to rather inconsistent results. The critical conductance was studied in |21]. We will 
discuss some of these results later together with our own results. 

2.2. Conductance 

The two-terminal conductance of a particular sample is calculated from 

g = Trth, (3) 

where t is the transmission matrix. Its elements to/3 determine the transmission 
amplitude from channel a into channel (3. Two semi-infinite ideal leads are attached to 
the disordered sample and periodic boundary conditions are applied in the transverse 
direction. The Fermi energy is fixed at Ep = 0.01 and spin-orbit scattering is assumed 
to be absent in the leads. 

We use the algorithm of Pendry et al for our numerical calculations. The size 
of the sample varied from L = 20 to L = 200. We concentrate on the analysis of the 
L-dependence of the mean conductance in the neighbourhood of the critical point at 
Wc ~ 5.80. Because of the absence of self-averaging of the conductance, we need to 
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Figure 1. The mean conductance {g{L)) as a function of the system size L for various 
values of disorder strength W . Within the interval of system sizes 30 < i < 200, 
{g{L)) is almost constant for disorder W — 5.84. This indicates that finite size effects 
are absent and irrelevant scaling terms negligible. Solid lines arc fits to |(SJ) with critical 
parameters gc = 1.42, Wc = 5.838, v = 2.80 and A = -0.23. 



analyze data for a large number of macroscopically identical samples which differ only 
by the respective microscopic realization of the disorder. For each disorder strength, 
5.71 < W < 5.96, and each system size, we collected a statistical ensemble of A^stat > 10^ 
samples and calculate the mean conductance (g) and variance vsng = {g"^) — (g)'^. Here, 
(. . .) means averaging over a statistical ensemble. Mean values (g) with the uncertainty 
5g = {varg/Nsta.tY^'^ are used in the scahng analysis. As a typical mean value {g) 1.4 
and vaig ~ 0.36, the relative uncertainty Sg/ (g) of our data is of order of 0.0015. Typical 
numerical data are shown in figure Q and |21 

2. 3. Critical wave function 

Electron wave functions at quantum critical points of disordered systems exhibit a 
peculiar spatial structure j2nil2Zl- For our symplectic model, the critical wave functions 
for systems sizes of up to L = 260 were obtained by direct diagonalization using a 
Lanczos-algorithm. The fractality of the spatial fluctuations of the modulus of the 
normalized eigenstates V's(^) was determined from the power-law scaling of the g-th 
moment defined via a 'box-probability' 

N{1) 

P{qA) = Y.^ E \i^E{r)\y ^ y^'^y (4) 

The generalized fractal dimensions D{q) = T{q)/{q — 1) or the so called /(a(g))-distri- 
bution was thus derived [211 1211 ED], where fli{l) is the the i-th box of size / = XL. The 
r(g) and f{a{q)) are related by a Legendre transform. 
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2.4- Scaling analysis 

At the critical point W = Wc, the conductance does not depend on the system 
size. In accordance with the scahng theory of locahzation j^I], we assume that in 
the neighbourhood of the critical point the L dependence of g is governed by the critical 
exponent z/, 

g^W, L)=g, + A{W - W,)L^'\ (5) 

We collected numerical data g{W,L) for more than = 200 values of W and L and 
obtain the critical parameters, gc-, Wc and v from a fit of the numerical data to the 
scaling ansatz (0). We made sure that more sophisticated fits [221 which include higher 
order terms in disorder are not needed because our data exhibit already a perfect linear 
ly-dependence (see figure E}, Also, subject to the uncertainty of our raw data, we found 
that possible irrelevant scaling fields IHH] are weak and do only marginally influence 
our results. 

The most simple analysis of critical parameters can be done as follows: first, we 
calculate parameters of the linear ly-dependence of the conductance 

g{W,L)=g^{L) + Wg^{L). (6) 

Comparing © with (0), we see that the slope gi depends on the system size as 

9i{L) oc L'l\ (7) 

Then, a power-law fit gi{L) vs. L gives us the critical exponent. 
We also minimize numerically the function 

F=ljT.(M-2MW,L))-gUW,L)f (8) 
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Figure 3. The scaling parameter A calculated for quasi-ld geometry versus disorder 
strength for varying system width M. The data confirm a critical disorder Wc ~ 5.835 
and Ac ~ 1.87 ±0.02. 



with respect to unknown parameters Wc, v and Qc and A. To estimate the accuracy of 
the final result, we repeat Nra times the minimization of the function F with input data 
{giyV^V)) randomly fiuctuating within their respective error bars. To understand the 
role of the finite size effects, we use in the scaling analysis only data for systems with 
sizes Lmin ^ L < Ljnax and study how the obtained critical parameters depend on the 
choice of and Lmax- 



3. Results and discussions 

3.1. Critical parameters 

Our analysis showed that neither the critical disorder nor the critical conductance 
depend on the choice of the minimal and maximal system size, Lmin and L^ax- From all 
the data we concluded that 

Wc = 5.838 ± 0.007 (9) 

The estimate of the critical disorder differs considerably from the result of Fastenrath et 
al 112] {Wc = 5.74), obtained from scaling analysis of the localization length of quasi-ld 
systems. Since previous results were derived from small system's widths, the discrepancy 
could be explained as an effect of finite size corrections. To support this assumption, 
we performed numerical simulations for quasi-ld systems of size M x L. Figure El shows 
the calculated data for the ratio A = \m/M, the localization length Am divided by the 
system width M jSHEll- For our purposes, we run A with an accuracy of only 1%. This 
was sufficient to confirm that our estimation of the critical disorder derived from the 
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Figure 4. The critical exponent v obtained from numerical data with L^^in < L < 
Lmax- Full symbols shows how v depends on the choice of imax for fixed minimal 
system size imin- Open symbols show how ly depends on Lmin for fixed Lmax, 
Lmax = 140 (O), and L^ax — 200 (A). The shaded area highlights our estimate 
of the critical exponent — 2.80 ± 0.04. 



conductance calculations indeed coincides with the one obtained from the locahzation 
length. Also, we get an estimate of the critical value of the scaling parameter A, 

Ac -1.87 ±0.02 (10) 

which agrees well with a recent result of Asada et al 1221 obtained for a SU(2) model. 

From our conductance data, the critical exponent v was obtained for the first time 
using finite size scaling of the calculated electrical two-terminal conductance, which is, 
unlike the localization length, an easy way to measure physical quantity. With the finite 
size scaling analysis described in section ()2.4|1 . our estimation of the critical exponent 
gives 

z/ = 2.8 ±0.04, (11) 

which is close to a numerical result v = 2.746 ± 0.009 for a SU(2) model published 
recently j^, but considerably larger than u = 2.05 ± 0.08 [H] and smaller than 
u = 2.88 ± 0.15 as obtained in j2l]. In contrast to the critical conductance and critical 
disorder, the critical exponent is sensitive to the system size as can be seen in figure 0] 
The estimate of u is therefore highly non-trivial and error bars obtained from various 
fit procedures differ. 



Critical conductance and wavefunction of two dimensional Ando model 

0.7, . ^ . ^ . ^ . — 





L= 82 


1.418 


< 


L=100 


1.415 


■ 


L=120 


1.417 


o 


L=140 


1.418 




L=200 


1.417 




W=5.84 



Figure 5. Critical conductance distribution p{g). Inset: the derivative dp{g)/dg 
which shows (/non-an to bc slightly larger than 2. 



3.2. Critical conductance 

The critical two-terminal conductance obtained by us for square samples with periodic 
boundary conditions in the transverse direction, 

= 1.42 ± 0.005, (12) 

can be compared with the result for a SU(2) model studied recently in [5^. From the 
corresponding numerical data for size L < 97 (T. Ohtsuki, private communication) we 
calculated lim^^^oo fi'(-^) = 1-411, which is in very good agreement with our present 
result. 

Assuming for square systems a = g/L^'^~'^\ we have at the critical point gc = 
As mentioned above, for L > 20 our critical conductance data do not exhibit any finite 
size effect. We do not expect that the numerical analysis causes any inaccuracy of the 
estimation of Qc. Still, when comparing with ctc, we have to keep in mind that small 
corrections due to the properties of leads (we assume perfect leads without spin-orbit 
scattering) might be responsible for a small difference between and gc- However, we 
do not expect this difference to be larger than ~ and neglect it for further purposes. 

Figure |31 shows the critical conductance distribution p{g). We present data for 
disorder W = 5.84 and system size 82 < L < 200. We see that indeed p{g) does not 
depend on the system size at the critical point. Besides the shape of the distribution, 
which is known also from previous studies, we want to point out the non-analytical 
behaviour of the distribution at the point g = (y'non-an which is slightly larger than 2. 
This agrees with a qualitative estimation of Muttalib et al although their work 
concerns quasi one dimensional systems with orthogonal symmetry. 




Figure 6. The /(a(g))-distributions of 10 critical eigenstates at £^ « and 
= 5.83 for q = 0.0, ±0.2, ±0.5, ±0.8, ±1.0, ±1.5, ±2.0, and ±3.0. The parabolic 
fit is determined by ao = 2.107 ± 5 • 10~^. The inset shows the scaling of lnP(l,A) 
vs. logarithm of the box size A. Data are obtained from ten L = 260 samples with 
different disorder realization. The dashed line has steepness t(1) — 1.889. 



3.3. Fractal dimensions and Fal'ko-Efetov relation 

In figure inithe /(Q;(g))-distributions of 10 critical eigenfunctions of systems with L = 260 
are sliown togetlier witli a parabolic fit f{a{q)) = d — («(g) — a;o)^/[4(ao — d)] with 
ao = (^{q = 0) = 2.107 ± 0.005. This value is lesser than the one obtained previously 
for smaller sample size L = 150 9J. The inset exhibits the logarithm of P{1, A) plotted 
versus In A for 2 < I < 25. The data belonging to critical eigenstates obtained from 10 
disorder realizations clearly obey power-law scaling. 

From the multi-fractal analysis of critical eigenfunctions of L = 260 samples, we 
obtained an information dimension D{1) = 1.889 ±0.001 using 10 different realizations 
with disorder strength W = 5.83. With this data, we checked the validity of a formula 
suggested by Fal'ko and Efetov [T7j which should hold in the limit 1 < gDh, 

where 13 = 1/2, 1 and 2 for or orthogonal, unitary and symplectic symmetry, 
respectively. With a = e^gD (diffusion constant D and density of states q) we get 

q 

27r/3[2-D(g)]l' 

Taking g = 1, we found our results for (Jc and -D(l) to be in very good agreement 
with ()14p . Therefore, this relation is fulfilled also for the symplectic symmetry in a 
similar manner as found previously for the quantum Hall case ISH]- In comparing 
these results one has to keep in mind that in the present paper two spin channels 
were considered, i.e., the conductance has to be divided by 2, whereas in the QHE case 
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Figure 7. The generalized fractal dimension D{q) of a critical eigenfunction as a 
function of q calculated via the box-probability Q. The linearity D{q) = 2 — q/k 
holds for q < 1.5, with k = 0.112 leading to = 0.71 e'^/h. 



only one spin direction was taken into account. Note also that the value of the disorder 
used in the multi-fractal analysis, W = 5.83, differs slightly form the critical disorder 
Wc = 5.838. Therefore, the conductance g{W = 5.83) ~ 1.43 mentioned in ref. jSH] 
is a little larger than Qc- The validity of the proposed linear relationship ()13|1 holds at 
least for g < 1.5 which can be seen in figure [3 where D{q) calculated according to (jl)) is 
plotted versus q. 



3.4- Comparison with Janssen's formula 

A formula that connects fractal properties of critical eigenstates in square samples with 
the finite size scaling variable A = Xm/M of quasi-ld systems was conjectured by 
Janssen fHlE]- At the critical point, Ac is scale independent and depends only on the 
Lipschitz-Holder exponent and the spatial dimension of the system d = 2 

Ac= , ^ ,y (15) 

TT[ao — d) 

Equation (fT3|) was reported to hold for 2d disordered systems in the quantum Hall 
regime fHllSZI- In the present symplectic system, however, it is not met by our results. 
Also, assuming the parabolic approximation to be valid, we can use the relation between 
D{1) and ao, 

D{1) = a{l) = A- ao, (16) 
and connect (jl4|) and (jl5|) to obtain a relation between A^ and the critical conductivity 
A 

Again, contrary to what is observed in the QHE regime, the relation (fTTj) is not satisfied 
for a symplectic system described by the Ando model. A possible reason for this failure 
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might be that in the present 2d model, the scahng parameter A,, may not represent the 
typical localization length of the system as has been assumed in |18j . 

4. Conclusion 

We studied the electrical two-terminal conductance and the spatial fluctuations of 
electron eigenf unctions near the metal-insulator transition of the two dimensional Ando 
model. Using finite-size scaling, we obtained for this symplectic model a critical 
exponent u = 2.8 ± 0.04, which governs the size dependence of the conductance in the 
critical regime, and a critical disorder Wc = 5.838 ± 0.007. Our results for the critical 
conductance and for the fractal dimension of critical eigenstates confirm the validity 
of Fal'ko and Efetov's prediction (fT^ . However, comparison of the fractal dimensions 
with the critical value of the finite size scaling parameter Ac, calculated for quasi one- 
dimensional systems, indicates that Janssen's formula (fTHjl is not fulfilled by our data 
in the present model. 

Unlike the localization length, the two-terminal conductance should be easily 
accessible in experiments. Our numerical confirmation of the relation between fractal 
dimensions of critical eigenf unctions and the critical conductance provides an additional 
argument for universality at the metal-insulator transition in 2d symplectic models. 
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